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Abstract 

The leading chiral contributions to the quark and gluon components of the 
proton spin are calculated using heavy-baryon chiral perturbation theory. 
Similar calculations are done for the moments of the generalized parton distri- 
butions relevant to the quark and gluon angular momentum densities. These 
results provide useful insight about the role of pions in the spin structure 
of the nucleon, and can serve as a guidance for extrapolating lattice QCD 
calculations at large quark masses to the chiral limit. 
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The spin structure of the proton has been an active subject of investigation in nuclear 
and high energy physics since the EMC (European Muon Collaboration) measurement of the 
gi{x, Q^) structure function in polarized muon deep-inelastic scattering EMC found that 
a surprisingly small fraction of the proton spin resides in quark spin — in flat contradiction 
with the well-known SU(6) quark model prediction 0. In the last decade, many follow- 
up experiments have been performed at the CERN, SLAC, and DESY laboratories, and 
the next-to-leading order QCD analyses of these data confirm the original discovery of 
EMC 0. The outstanding question is then what are the other contributions to the proton 
spin? According to quantum chromodynamics (QCD), the angular momentum of a strong- 
interacting system can be written as a sum of quark spin (S^), quark orbital (Lg), and gluon 
angular momentum (Jg) contributions [Q, 

Jqcd = / c^^x^V^St/' + J £x^^x X (-iD) 

+ J c/^xx X (E X B) , (1) 

where D = 9 + igA is the covariant derivative, A the gluon potentials, ijj the quark fields, 
and E and B the chromoelectric and magnetic fields. Summations over flavors and colors are 
implicit. Using the above operator, one can write down a gauge-invariant, frame-independent 
decomposition of the proton spin 0] , 

I = lAIlif,') + L,{^,') + J,{^^') , (2) 

in a class of frames in which the proton has a definite helicity 0. Individual terms in 
the above equation are defined as the expectation values of the corresponding operators in 
the proton helicity eigenstate. Because these operators are not separately conserved, their 
matrix elements depend on renormalization schemes and the related scale /i^. 
Recent analyses of the polarized deep-inelastic scattering data yield f^, 

AE(/i2 ~ SGeV^) = 0.16 ±0.08 . (3) 

Meanwhile, it has been realized that the total quark and gluon contributions to the proton 
spin, Jq = AS/2 + Lq and Jg, can be obtained from new sum rules involving generalized 
(off-forward) parton distributions |Q. These new distributions are found to be accessible 
experimentally in deeply virtual Compton scattering or exclusive meson production [Q. 
On the other hand, the matrix elements AS, Lg and Jg can be calculated theoretically in 
lattice QCD and QCD sum rule methods 0. In the former case, the up and down quark 
masses used in the present simulations are much larger than the physical ones, and the 
corresponding pion mass is larger than 500 MeV. Therefore, it is crucial to make reliable 
extrapolations from the lattice data down to the physical pion mass. 

In this paper, we study the pionic contribution to the spin structure of the proton. This 
contribution can be calculated because pions are Goldstone bosons from the spontaneous 
breaking of the chiral symmetry present in the QCD lagrangian. The explicit breaking 
of the symmetry through small up and down quark masses induces small physical pion 
masses. Effective field theory technology allows study of the symmetry breaking effects in a 
systematic expansion of m7r/47r/^. Here we focus on the leading nonanalytical contributions 
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of the type m^lnm^ to the quark and gluon components of the proton spin, Jq^g. The total 
nonanalytic contribution to the proton spin, however, must vanish because the total proton 
spin is independent of the pion mass. 

What we actually calculate below is the chiral contribution to the proton form factors 
of the quark and gluon parts of the QCD energy momentum tensor T^^. A simple counting 
shows that they have three form factors H,|H], 



P 



P 



u{p' 



A 



1,9 



(4) 



where p = {p+p')/2, q = p'—p, and (• ■ ■) means that the indices enclosed are made symmetric 
and traceless. Aq(g)(0, /i^) is the momentum fraction of the proton carried by quarks (gluons), 
and therefore ^^(O,;^^) + Ag{0, fi"^) = 1. The quark (gluon) angular momentum contribution 
to the proton spin is M 



(5) 

0. An explicit 



Thus the sum rule Jg(/i^) + JgifJ''^) = 1/2 implies Bq{0,fi^) + _Bg(0,/i^) 
one-loop verification of this general result in QED can be found in Ref. 0. 

Since we are interested in the form factors A{q'^,fi^) and B{q'^,fi^) at zero momentum 
transfer, chiral perturbation theory is a legitimate tool to use. As a low-energy effective field 
theory of QCD, chiral perturbation theory treats chiral symmetry and the symmetry break- 
ing patterns of QCD in the most general fashion. In the single nucleon systems, the baryon 
masses are considered heavy, and the heavy-baryon chiral perturbation theory (HBxPT) 
has been formulated to restore the systematic power counting [|10|. In this formalism, we 
rewrite the form factor in the Breit frame as 



P 



P 




q^ 



V 4M2 
+ {Aq,giq',^^') + Bq,giq',^^'))v^>^ 
+ C,,,(g^/i2)gWM}A^(^;) , 



4M2 



S''\S 



(6) 



where the proton velocity v = {p + p')/2M, polarization vector 5"^ = i'y^cTf^i^v'^ reduces to 
(0,(t/2) in the = (1,0,0,0) frame (A^A^ = 2M). In the following, we will focus on the 
spin-dependent term (the second in the bracket), and compute the chiral corrections to its 
coefficient, Ag^g{0) + S,,g(0) = 2Jg,g. 

According to the formalism developed in Refs. 



11|, the quark and gluon energy- 



momentum tensor operators are "matched" onto a sum of hadronic operators with the 
identical quantum numbers, including pure pionic operators, single baryon operators as well 
as multiple baryon operators. 



1 rp(fiv) 

+ ^ ('?,9)A 



(/i',A^) + 



(7) 



where A^ is a scale at which HBxPT is applied. Each of these operators contains an infinite 
number of terms organized according to the HB^PT power counting. For instance, the 
leading terms in the pure pionic operator are 
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FIG. 1. Feynman diagrams contributing to the leading chiral behavior of the fractions of the 
nucleon spin carried in quarks and gluons. The dashed hnes represent pions. The diagram (d) 
denotes the wave function renormahzation contribution. 

tIZW^K) = afe,).(/i^)f Tr[a(^St9^)S] + 6(,,,).(/.2)|a(^Tr[St5^)S + S5'^)St] + ■ ■ ■ 

(8) 

where S = -u^ = exp[i7r°r"//7r], n"" is the pion field with isospin index a, the isospin 
Pauh matrices, /^r the pion decay constant = 93 MeV. The physical meaning of a{q,g)n is the 
momentum fraction of the pion carried by quarks (gluons) in the chiral limit. Hence, 

aqnifJ'^) + %7r(/i^) = 1 • (9) 

Experimentally, it is found that 

a,^(/i2 ~ 1 GeV^) = {x,)^ - 0.5 (10) 



for a physical pion [|T2| . The leading pionic operator contributes to the leading chiral behavior 
of the spin fractions Jg g shown in Fig. 1(c). 

The leading terms in the single nucleon operator can be written as 

+C(,,,)^(/^2)iV^(/^S'^)iV + ... (11) 

where = i {v)d'^u — ud^u^^ = —FT^r'^d'^ir'^ + ■ ■ ■. The second term contributes to Jq^g 
through tree and one loop diagrams (Fig. 1, diagrams (a), (b), and (d)). In fact, b(q^g)N/2 
is just Jq^g in the limit where quark masses vanish 

Jq,gif^') = &fe.)^(/i')/2 , (12) 
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(a) (b) (c) 



FIG. 2. Same as Fig. 1. Contributions from the delta intermediate states (double lines). The 
diagram (c) denotes the wave function renormalization contribution. 

where the superscript indicates the chiral hmit. The total nucleon momentum and spin in 
the chiral limit impose the constraints 

agAf(/^^) + agAf(/i^) = 1 , 

+ = 1 ■ (13) 

Combining the result from different diagrams in Fig. 1, we find the following leading 
chiral logarithm contributions to the spin structure of the nucleon, 

^.,.(/^') = \ { + 3[afe,).(/i^) - h.,9)N{^^')]^^,^^ (^) } ^ ■ ■ ■ ' ^^^^ 

This is the main result of this paper. The chiral logarithms in Jq{i^'^) + Jgifi''^) cancel because 
of the constraints in Eqs. |TB|). 

The delta resonance plays a special role in the nucleon structure physics. In a world 
where the number of quark colors (A^^^) go to infinity, the delta resonance is degenerate with 
the nucleon, i.e. A = Ma — Mn = 0{1/Nc). If, at the same time, the standard spontaneous 
chiral symmetry breaking occurs, the chiral corrections to the nucleon properties are strongly 
affected by the degenerate delta contribution. In the real world where A ~ 300 MeV, the 
significance of the delta contribution depends strongly on the spin-isospin channel to which 
an observable belongs |TB[ . It can be calculated in an effective field theory approach in which 



A is counted as the same order as 10,1 



In Fig. 2, we have shown three different possible ways that the delta resonance can 
contribute as an intermediate state to the spin structure of the nucleon: First, the delta is 
present when the energy-momentum tensor of the pion contributes (Fig. 2b); second, the 
delta contributes to the wave function renormalization of the nucleon (Fig. 2c); and finally 
the energy-momentum tensor of the delta contributes directly (Fig. 2a) 

= -3h,,o)Aif^'W {^'v^' [S''\ S^] A,} + ... (15) 

Here the leading term contributes to the delta spin in the chiral limit. Indeed, the total spin 
of the delta leads to the following constraint. 



&gA(/i^) + &9A(/i^) = 1 



(16) 



The calculation of the delta contribution to Jq^g^jJ^) is straightforward. The combined 
leading non-analytic chiral correction is 



1 
2 



(1 7D 



9j /2\ o /2\ 



where 

K{m^,A) = {ml -2A'') log 



J (3 -471^)2 



A 



+ 2AJA2 -m2 log 



A^ — ml + 



(17) 



A + W A2 - m2 + ie 



(18) 



Again, all the delta contributions cancel when summed over quark and gluon operators 
because of the constraints on the low-energy constants. 

Another nontrivial check of the above result is to let the number of colors Nc go to infinity. 
In this limit, the individual nucleon and delta contributions scale as Nc- However, because 
2^/2g^N/^/3 = qa and bi^q^g)^ = &(g,g)A at large Nc [|l5l, the order Nc chiral contributions 
cancel. In the real world, we can use the large- iVc approximation of the delta observables to 
get an estimate of the total chiral contributions. 



.2^ _ TO , ^9 



JqAf^') = j'q + 



(47r/.) 



q/n 



K{mT,, A) — ml In 



X 



(19) 



which depends on the unknown J°. It is interesting to note that the non- analytic corrections 
will be largely cancelled if J° ~ {^)q/TT /2, a possibility favored by empirical considerations. 
In this linear ml extrapolation of lattice calculations to the chiral limit might be 

justified. 

In Refs. [|l^], the quark and gluon angular momentum distributions are found to be 
Jq{x) = x(J2{x) + Eq{x, 0, 0))/2 and Jg{x) = x{g{x)+Eg{x, 0, 0))/2, where S(a;) and g{x) are 
unpolarized (singlet) quark and gluon distributions and Eq^g{x, 0,0) are generalized parton 
distributions . The moments of these distributions are related to the form factors of the 
twist-two operators P]. In leading order in q, 



^p'|0WM2...Mn|p^ = Ar,fl{0, fi^)N{v)v^''v^\..v^'"N{v)M' 



n-l 



+ [An,oiO, /i^) + 5„,o(0, fi^)]N{v)[S''\S ■ q]v^'\..v>'-N{v)M 



n-2 



+ ... 



Then the relation is 



/■I 1 

j_^x''-^Jq,g{x,fx^)dx = -[A„,o(0,/i^) + 5„,o(0,/x^)] , 
for n = 2,4, .... By going through the same calculation as above, one can show that. 



(20) 



(21) 



X" Jq,g{x,IJ, )dx = bn{q,g)N{lJ' ) 1 



3-M.ln 



- -b 



'n{q,g)N 



15, 



,s)A(/i )j 



(47r/.) 



(3 • 4nU 



-K{m^,A) + -- 



(22) 
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for n = 4, 6, where bn{q,g)N and hn^q^g-^/^ are the coefficients of the nucleon and delta spin- 
dependent twist two operators in the effective theory. Apart from the absence of the pion 
contribution, the other terms are similar to the n = 2 moment. 

To summarize, we have calculated the leading chiral contribution to the spin structure 
of the nucleon. These results provide interesting insight about the role of the pion in the 
composition of the nucleon spin. It also provides useful guidance for extrapolating lattice 
QCD calculations at large quark masses to the chiral limit [0,0. 

This work is supported in part by the U.S. Dept. of Energy under grant No. DE-FG02- 
93ER-40762-245. 
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